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One specific hypothesis explaining the evolution of extra-pair reproduction (EPR) by socially monog- 
amous females is that EPR is under indirect selection because extra-pair offspring (EPO) sired by 
extra-pair males have higher additive genetic value for fitness than the within-pair offspring (WPO) a 
female would have produced had she solely mated with her socially paired male. This hypothesis has 
not been explicitly tested by comparing additive genetic value between EPO and the WPO they replaced. 
We show that the difference in additive genetic breeding value (BV) between EPO and the WPO they 
replaced is proportional to the genetic covariance between offspring fitness and male net paternity gain 
through EPR, and estimate this covariance with respect to offspring recruitment in free-living song spar- 
rows (Melospiza melodid) . Recruitment and net paternity gain showed non-zero additive genetic variance 
and heritability, and negative genetic covariance. Opposite to prediction, EPO therefore had lower BV for 
recruitment than the WPO they replaced. We thereby demonstrate an explicit quantitative genetic 
approach to testing the hypothesis that EPR allows polyandrous females to increase offspring additive 
genetic value, and suggest that there may be weak indirect selection against female EPR through reduced 
additive genetic value for recruitment of EPO versus WPO in song sparrows. 
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1. INTRODUCTION 

The forces driving the evolution of extra-pair reproduction 
(EPR) by socially monogamous females, and of polyandry 
more generally, remain widely debated and unresolved 
[1-8]. One major hypothesis is that female EPR is under 
indirect selection because the extra-pair offspring (EPO) 
produced through EPR are fitter than the within-pair 
offspring (WPO) a female would have produced had she 
solely mated with her socially paired male [1-3]. This 
hypothesis is often implicitly genetic; in many socially mon- 
ogamous but genetically polyandrous species, extra-pair 
(EP) sires do not obviously provide additional resources 
to a female's EPO. Any intrinsic difference in fitness 
between a female's EPO and WPO that could underlie 
indirect selection for female EPR is therefore postulated 
to reflect increased additive and/or non-additive genetic 
value of EPO, and hence increased value of the EP sire 
compared to the female's socially paired male with respect 
to one or both components of genetic variation (commonly 
termed 'good genes' and 'compatible genes' theories, 
respectively [1,2]). Most explicitly, if females produced 
EPO of higher additive genetic value for fitness than the 
WPO they replaced (reflecting higher additive genetic 
value of the EP versus within-pair (WP) sire) and there 
was additive genetic variation in female propensity for 
EPR, genes underlying female EPR could become 
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correlated (in linkage disequilibrium) with genes under- 
lying high paternal or offspring fitness. Female EPR 
could then evolve through indirect selection [3,9] . 

Such genetic covariances are difficult to measure, 
especially in free-living populations. One widely implemen- 
ted short-cut to testing the broad hypothesis that female 
EPR is under indirect selection is therefore to compare 
the phenotypes of maternal half-sib EPO and WPO from 
the same clutch, brood or litter [2-4]. Such comparisons 
have the advantage that phenotypic effects of natal environ- 
ment and maternal genes may be relatively similar across 
maternal half-sibs. Overall effects of paternity (and hence 
paternal genes) on offspring phenotype can therefore be 
estimated independent of environmental and maternal 
variation that can severely confound phenotypic com- 
parisons among sires or broods [3,10]. Consequently, 
numerous field studies have compared phenotypic traits 
such as body size, immune response and early survival 
between maternal half-sib EPO and WPO from the same 
brood, and recent studies have compared major fitness 
components, including lifespan and lifetime breeding 
success [2,4,8,10-12]. Such comparisons are valuable in 
quantifying overall phenotypic differences between EPO 
and WPO, but face certain limitations in testing the specific 
hypothesis that females increase the genetic value of their 
offspring through EPR. 

Explicit test of this hypothesis clearly requires the 
genetic values of WPO and EPO, rather than solely 
their phenotypic values, to be quantified and compared. 
Comparing phenotype, even between maternal half-sibs, 
could provide misleading conclusions regarding relative 
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genetic value if environmental or maternal effects on indi- 
vidual phenotype covary with paternity within a brood or 
litter. For example, EPO may be conceived early within a 
brood, perhaps reflecting the timing of EP copulations 
[8] . Early offspring may also have higher fitness than later 
offspring independent of their paternity, owing to variation 
in competitive ability resulting from hatch or birth order 
within a brood. EPO may therefore be fitter than WPO 
because of their brood position rather than because they 
are of intrinsically higher genetic value [8, 1 3, 1 4] . The con- 
verse scenario, where EPO hatch later and have low fitness 
owing to their brood position, could also be postulated. 
Indeed, phenotypic comparisons between maternal half- 
sibs from the same brood do not control for any form of 
within-brood variation in environmental or non-genetic 
parental effects on offspring fitness in relation to paternity, 
whether owing to brood position or differential parental 
investment or any other mechanism [2,3,8,10,14]. Such 
effects could cause the difference in genetic value between 
EPO and WPO to be overestimated based on compari- 
sons of phenotype, or underestimated if environmental 
effects systematically counteract genetic effects, providing 
misleading evidence regarding the benefit of EPR to 
polyandrous females. 

This situation reflects a general limitation of testing 
the hypothesis that female EPR is under indirect selection 
by comparing phenotypes between observed maternal 
half-sib EPO and WPO. Strictly, the indirect selection 
hypothesis predicts that EPO will be fitter than the 
WPO a female would have produced had she solely 
mated with her socially paired mate (i.e. the WPO that 
the EPO replaced). It does not strictly predict that EPO 
will be fitter than the WPO a female produced anyway, 
despite having produced the EPO. However, the pheno- 
types of WPO that were not produced clearly cannot be 
measured and compared with the phenotypes of the 
EPO that were produced instead. The approach of com- 
paring phenotypes between observed maternal half-sib 
WPO and EPO therefore requires an additional assump- 
tion: that the WPO that were produced are similar to the 
WPO that would have been produced if the female had 
not produced EPO instead. This may be valid (on aver- 
age) with respect to genetic value, but it may not be 
valid with respect to environmental effects on phenotype, 
whether reflecting brood position or parental investment 
or any associated maternal effects. 

A further constraint on comparing observed phenotypes 
between maternal half-sib EPO and WPO from the same 
brood or litter is that such analyses are necessarily restricted 
to mixed-paternity broods and cannot include broods that 
comprise entirely EPO (where EPR occurred to the degree 
that the female's socially paired male was completely cuck- 
olded [10-12]). In these cases, there are no WPO whose 
phenotypes can be measured and compared with those of 
the EPO. Broods where all offspring were sired by an EP 
male, which may be important drivers of selection on 
EPR and highly informative in understanding the basis of 
such selection, are therefore excluded from key analyses. 

Given the specific aim of testing the hypothesis that 
polyandrous females increase the genetic value of their 
offspring through EPR, these problems of distinguishing 
genetic and environmental effects underlying the relative 
fitness of WPO and EPO, comparing EPO with the 
WPO they replaced and therefore were not produced 



and incorporating broods that comprise entirely EPO 
can, in theory, be resolved by explicitly estimating quantita- 
tive genetic rather than phenotypic differences between 
maternal half-sib EPO and WPO . For example, an individ- 
ual's breeding value (BV) for any trait, defined as the sum of 
the average additive effects of its alleles on that trait (hence 
measuring additive genetic value), can be estimated given 
sufficient phenotypic data from known relatives spread 
over environments and parents [15,16]. This approach 
allows BVs to be estimated for individuals whose own 
phenotype was not observed, as commonly done for 
sex-limited traits in individuals of the non-expressing sex. 
The same approach can, given sufficient observations of 
relatives' phenotypes, allow estimation of BVs for individ- 
uals that do not exist, such as the hypothetical WPO a 
female would have produced had she not produced EPO 
instead. The differences in BVs and phenotypic values 
between maternal half-sib WPO and EPO can then be 
compared, informing the degree to which differences 
in additive genetic value can be inferred from standard 
phenotypic comparisons. 

Although quantitative genetic analysis potentially pro- 
vides the most explicit test of the hypothesis that a 
female's EPO are of higher additive genetic value than the 
WPO they replaced, this approach imposes formidable 
data requirements, is vulnerable to bias and has not been 
implemented in a wild socially monogamous but geneti- 
cally polyandrous population. We derived a general 
quantitative genetic expression describing the difference 
in BV between a female's EPO and the hypothetical 
WPO they replaced. We then used comprehensive genetic 
pedigree data from free-living song sparrows (Melospiza 
melodid) to parametrize this expression and explicitly test 
the hypothesis that a female's EPO have higher BV for 
survival to recruitment, one major fitness component, 
than the hypothetical WPO they replaced. We compared 
the estimated differences in sex-specific BVs and phenoty- 
pic values for recruitment between maternal half-sib EPO 
and WPO and consider the degree to which conclusions 
are consistent. We thereby formulate and implement an 
approach to explicitly comparing additive genetic value 
between EPO and the WPO they replaced, and discuss 
the contribution of additive genetic variation in offspring 
recruitment to indirect selection on female EPR. 



2. MATERIAL AND METHODS 

(a) Estimating the difference in breeding value 

'Animal models' are generalized linear mixed models where 
pairwise coefficients of kinship (k) calculated from pedigree 
data define a matrix proportional to the variance -covariance 
matrix of additive genetic random effects [15,16]. They allow 
additive genetic variances {V^) and covariances (cova) in and 
among phenotypic traits to be estimated given unbalanced 
data, inbreeding and selection [15,16]. They also allow BVs 
for focal traits to be estimated for all individuals included in 
the pedigree [15,16]. BVs for hypothetical individuals that did 
not exist (such as the WPO that EPO replaced) could therefore 
be directly estimated by coding these individuals in the pedigree, 
or calculated as half the sum of their parents' estimated BVs. 
The mean difference in BV between a female's EPO and the 
WPO they replaced (£[ABV]) could then be calculated. 

However, estimated BVs may be biased if appropriate 
environmental effects on mean phenotype, environmental 
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variances and genetic groups are not modelled, and error 
and uncertainty may be substantial if phenotypes are observed 
for only few relatives of each individual [17-19]. An individ- 
ual's estimated BV can then largely reflect its own observed 
phenotype and consequently be biased by environmental or 
non-additive genetic sources of phenotypic variation. Further- 
more, post hoc comparisons of estimated BVs among groups 
of individuals are likely to be biased if the groups are not 
specified in the underlying animal model, and associated stat- 
istical confidence is hard to estimate [18]. The most rigorous 
and unbiased approach is therefore to estimate i?[ABV] 
directly within an animal model rather than by post hoc 
comparison of estimated BVs of EPO and WPO [18]. 

Direct estimation of i?[ABV] requires this difference to be 
expressed as a genetic covariance between phenotypic traits 
that can be explicitly specified within an animal model. Alge- 
braic transformation (electronic supplementary material) 
shows that £[ABV] = i(cov A (N E - N C ,W)/E[N B ]), where 
W is fitness, N E is the number of EPO a male sired and 
N c is the number of offspring that male lost through cuckol- 
dry (the number of EPO produced by the male's socially 
paired female). The term (AT E — N c ), therefore, describes a 
male's net paternity gain (or loss) through EPR, and E[N E ] 
is mean N B across all males. i?[ABV] can therefore be 
estimated as half the genetic covariance between (iV E — Nc) 
and W across all males, divided by E[N E ] . E[N E ] can be 
calculated directly from paternity data and the genetic covari- 
ance cov A (iV E — N C ,W) can be estimated from a bivariate 
animal model for (N E — N c ) and W. The mean difference in 
BV between a female's EPO and the WPO they replaced can 
thus be estimated from a consistent model without post hoc 
analysis of estimated BVs. Positive i?[ABV| would indicate 
that females' EPO had a higher BV for fitness than the 
hypothetical WPO they replaced, on average. 

This basic expression for .E[ABV] assumes that additive 
genetic effects on fitness are the same in males and females 
and hence that the inter-sex genetic correlation for W is 
r mf ~ 1 . It also assumes that the genetic correlation between 
a male's net paternity gain through sons and daughters is 
ca 1 . The latter constraint can be relaxed by considering net 
paternity gain through sons and daughters as separate traits. 
The mean difference in BV between a female's EP sons and 
the WP sons they are assumed to have replaced can be esti- 
mated as £[ABV S ] =i(cov A (N ES - N CS ,W)/E[N BS ]), where 
iV ES and N C s are the numbers of sons that a male gains and 
loses through EPR. An equivalent expression holds for daugh- 
ters, such that £[ABV D ] = ±(cov A (iV ED - N CV ,W)/E[N EU ]) 
(electronic supplementary material) . 

(b) Study system 

Direct estimation of i?[ABV| requires AT E , N c and some 
measure of fitness W to be measured across numerous 
known relatives in a socially monogamous population with 
EPR. A long-term study of song sparrows on Mandarte 
Island, BC, Canada, fulfils these requirements. 

Song sparrows of both sexes typically breed two to three 
times per season starting from age one year, and are primarily 
socially monogamous with well-defined social pairings [20] . 
Females incubate clutches (typically three to four eggs) and 
both socially paired parents provision hatched offspring 
[20] . However, they are genetically polygynandrous, with fre- 
quent EPR [21]. Mandarte {ca six hectares) holds a resident 
song sparrow population, recently numbering ca 15-45 
breeding pairs, which has been studied intensively since 



1975 [20]. Each year, all breeding attempts are closely mon- 
itored and all nests located. All chicks surviving to ca 6 days 
post-hatch are marked with unique combinations of coloured 
rings to allow individual identification. The occasional immi- 
grants (1.1 yr _1 on average) are individually colour-ringed 
soon after arrival. This low immigration rate is sufficient 
to maintain allelic diversity [22] and prevent inbreeding 
from accumulating. All social parents of all offspring (those 
incubating clutches or provisioning chicks) are identified 
by observation. Owing to Mandarte's small size and the 
intensive fieldwork, the probability of resighting a surviving 
adult during any breeding season is effectively one. The 
local survival of all offspring from ringing at ca 6 days post- 
hatch to adulthood at age one year (hereafter 'recruitment') 
can therefore be documented with high confidence. The 
high local recruitment rate and absence of ringed sparrows 
breeding on surrounding islands suggest that emigration is 
relatively rare [12,20]. 

During 1993-2008, 99.4 per cent of ringed chicks and 
adults were blood sampled and genotyped at 13 polymorphic 
microsatellite loci to allow assignment of genetic parents 
[21]. Bayesian full probability models that incorporated gen- 
etic and spatial information allowed genetic sires to be 
assigned to 99.6 per cent of sampled offspring with at least 
95 per cent individual-level confidence, and revealed that 
ca 28 per cent of hatched chicks were sired by a male other 
than a female's socially paired mate and hence were EPO 
(compared with 24 per cent in a nearby mainland song spar- 
row population [21]). The probability of excluding a female's 
social mate as sire was 0.9998. Approximately 44 per cent of 
broods contained at least 1 EPO at ringing. All genetic 
mothers matched those assigned by social behaviour [21]. 
Sampled chicks were sexed using standard molecular 
methods [23]. Molecular sexes matched those attributed 
from reproductive behaviour for all recruited individuals. 
Male net paternity gain (N E — N c ) through all offspring, 
and through sons and daughters separately, was calculated 
for each adult (> 1 year old) male alive in each year during 
1993-2008. 

(c) Additive genetic (co)variances 

The expression E[ABV] = i(cov A (iV E - N C ,W)/E[N E ]) shows 
that non-zero ii[ABV] requires non-zero cov A (N E — N C ,W). 
This, in turn, requires non-zero V A in both (iV E — N c ) and 
W (currently measured as offspring survival to recruitment), 
because V A = 0 implies cov A = 0. Verbally, the hypothesis 
that EPO and WPO differ in BV for recruitment requires 
that both recruitment and male net paternity gain through 
EPR show non-zero V A and are heritable (h 2 > 0 [1,9]). 
Therefore, to estimate V A and h 2 in these traits, we first fitted 
univariate animal models to phenotypic data describing 
whether or not each offspring ringed on Mandarte during 
1993-2008 survived to recruit, and net paternity gain of all 
adult males alive in these years. 

Recruitment models included fixed effects of an off- 
spring's natal year and sex (because mean recruitment 
differs among cohorts and between males and females 
[20]), and a fixed regression on inbreeding coefficient 
(/, because recruitment shows inbreeding depression and 
failing to control for inbreeding depression can inflate esti- 
mates of V A [24,25]). Models also included fixed effects of 
EP status (EPO versus WPO) to minimize the degree to 
which any environmental variation associated with status 
could confound estimated genetic effects [17,18] and a 



Proc. R. Soc. B (2012) 



Genetic value of extra-pair offspring J. M. Reid & R. J. Sardell 1703 



status by sex interaction (because relative recruitment differs 
between male and female EPO versus WPO in mixed-paternity 
broods [12]). However, estimates of V A and h 2 remained quan- 
titatively similar when the latter two effects were excluded 
(electronic supplementary material). Random effects of an 
individual's natal brood, natal territory, mother identity and/ 
or social father identity were initially fitted to estimate var- 
iances owing to consistent effects of brood, location and 
parents providing care [16,19]. However, these variances 
were close to zero and estimates of fixed effects, V A and h 2 , 
remained quantitatively similar when they were excluded (elec- 
tronic supplementary material) . A bivariate animal model was 
then fitted to sex-specific recruitment data to estimate sex- 
specific V A , h 2 and inbreeding depression in recruitment and 
the inter-sex genetic covariance (cov mf ) and correlation (r mf ), 
given sex-specific fixed effects of year, /and EP status. 

Net paternity gain models included fixed effects of male 
breeding year, fixed regression on / and random effects of 
individual males to estimate 'permanent individual' variance 
(V m ) and account for repeated observations of individuals 
across multiple years. Random effects of a male's mother 
and social father identities were initially fitted, but estimated 
variances were close to zero and the estimates of inbreeding 
depression, V A and h 2 remained similar when these effects 
were excluded (electronic supplementary material). 

Bivariate animal models for offspring recruitment and male 
net paternity gain were then fitted to estimate cov A (iV E — 
N C ,W), with (iV E — N c ) measured across all offspring, and 
across sons and daughters separately. Appropriate fixed 
and random effects were fitted as above. 

(d) Analysts implementation 

Pedigree data for all offspring ringed during 1993-2008 
were compiled from genetic parentage data, taking the 
most probable sire for each offspring [21,26,27]. Pedigree 
data for 1975-1992 were compiled from observed social 
parentage and combined with the 1993-2008 genetic data 
to provide a full pedigree for all individuals alive during 
1975-2008. Assuming the unobserved EPR rate during 
1975-1992 was similar to the 28 per cent observed 
during 1993-2008, and that all mothers were correctly 
assigned by social behaviour, ca 86 per cent of all 1975- 

1992 pedigree links will be correct. Using the 1975-1992 
social pedigree data, even though uncorrected for EPR, is 
therefore likely to provide a more powerful analysis than 
simply using the corrected genetic pedigree data for 
1993-2008 (thereby assuming zero relatedness among the 

1993 breeders [26,27]). Inbreeding coefficients were calcu- 
lated from the 1975-2008 pedigree using standard 
algorithms [24,26]. Kinship between immigrants and exist- 
ing natives, and hence the / of offspring of immigrant- 
native pairings, was assumed to be zero relative to the base- 
line pedigree generation [24,28]. Inbreeding coefficients of 
immigrants are undefined relative to this baseline. Net 
paternity gain data from five immigrant males (with unde- 
fined /) were therefore excluded from analyses (electronic 
supplementary material) . 

Animal models were fitted using Bayesian methods 
implemented using MCMCglmm 2.14 in R v. 2.12.2 
[29,30]. Recruitment was modelled as a binary trait with 
binomial errors, logit link and fixed residual variance Pr = 1 
(by convention). Male net paternity gain was modelled assum- 
ing Gaussian errors (electronic supplementary material). Priors 
on fixed effects were normally distributed with mean zero and 



large variance (10 ). Parameter expanded priors on variance 
components used normally distributed working parameter 
priors with mean zero and variance 1000 and inverse-Wishart 
distributed location effect priors with degree of belief and 
limit variance of one (forming a scaled non-central F-distri- 
bution [15,30]). Prior genetic covariances were specified as 
zero. Posterior distributions were robust to reasonable variation 
in prior specifications (electronic supplementary material). 
Residual covariances between male and female recruitment, 
and between recruitment and male net paternity gain, were 
fixed to zero. Analyses used ^ 2 005 000 iterations, burn-in 
> 5000 and thinning interval ^ 2000 to ensure low autocorre- 
lation among thinned samples (less than 0.05). Posterior modes 
and 95% credible intervals (95% CI) for fixed effects, Va, h 2 , 
cov mf3 r mf an d cov A (iV E — N C ,W), were estimated across 
samples. The latent-scale heritability of recruitment was esti- 
mated as /zf atent = V A /(V A + 1 + 7^/3) given V R = 1 and 
logistic variance proportional to 77^/3 [31]. Aj 2 atent is interpret- 
able as the genetic intra-class correlation, or as the 
heritability of a latent variable describing an individual's liab- 
ility to recruit (where continuous variation is assumed to 
translate into recruitment at some threshold [15,31]). The 
probability-scale heritability of recruitment was estimated 
^ Ap r ob = (Va ■ X 2 /(l + /%) 2 )/((F A + 1) ■ X 2 /(l + Mr ) 2 + 
X ■ (1 — X)), where /j, r is the mean recruitment rate and X 
= /j,r/(1 + /ir) [31]. The heritability of male net paternity 
gain was estimated as A 2 at = V A /{V A + Vpi + Vr). 
The inter-sex genetic correlation for recruitment was 
estimated as r m f = coVmf /^/(J'Am x Vm)- The posterior 
distribution of the mean difference in BV between a female's 
EPO and the WPO they replaced (£[ABV]) was calculated 
by dividing the posterior distribution for cov A (AT E — N C ,W) 
by 2£[iV E ]. 

(e) Phenotypic variation 

Our secondary aim was to quantify the degree to which the 
sex-specific difference in BV between EPO and the WPO 
they replaced would have been accurately estimated by com- 
paring sex-specific phenotypic recruitment between EPO 
and WPO. However, the estimated difference in BV between 
EPO and the hypothetical WPO they replaced cannot be 
directly compared with the phenotypic difference between 
these individuals because hypothetical WPO have no obser- 
vable phenotype. To estimate relevant phenotypic variation 
as well as feasible, we fitted a linear model with binomial 
errors and the same fixed effects as the univariate animal 
model to phenotypic recruitment data for all offspring 
hatched in broods that contained at least 1 EPO during 
1993-2008. 



3. RESULTS 

(a) Recruitment 

Survival from ringing to recruitment was observed for 
2196 known-sex offspring hatched during 1993-2008, 
representing 807 broods (modal brood size 3, mean 
2.72, range 1-4). These broods included 234 known 
mixed-paternity broods and 124 EPO-only broods at 
ringing. The 2196 offspring comprised 796 WP daugh- 
ters, 312 EP daughters, 756 WP sons and 332 EP sons 
(totalling 1552 WPO and 644 EPO, and 1108 daughters 
and 1088 sons). 

Overall, 413 (18.8%) of the 2196 offspring recruited on 
Mandarte. The pedigree, pruned to the 2196 offspring and 



Proc. R. Soc. B (2012) 



1704 J. M. Reid & R. J. Sardell Genetic value of extra-pair offspring 



Table 1. Posterior modes (and 95% credible intervals) for additive genetic variance, heritabilities and inbreeding depression 
in survival to recruitment, and fixed effects of sex and extra-pair (EP) status, estimated across 2196 known-sex song sparrow 
offspring in a univariate animal model. h^ [ob was estimated assuming yii R = 0.19. Models without the EP status terms are 
provided as electronic supplementary material. 
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Table 2. Posterior modes (and 95% credible intervals) for sex-specific additive genetic (co)variances, latent-scale heritabilities 
and inbreeding depression in survival to recruitment, the inter-sex genetic correlation and fixed effects of EP status, 
estimated across 1088 male and 1108 female song sparrow offspring in a bivariate animal model. The large 95% CIs 
compared with table 1 reflect reduced sample sizes and power when males and females are considered separately. 
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Table 3. Posterior modes (and 95% credible intervals) for variance components, heritability and inbreeding depression in 
male net paternity gain (N B — N c ) estimated across 293 adult male song sparrows in a univariate animal model. 



additive genetic permanent individual inbreeding 

variance variance (V F {) residual variance (Fr) heritability (fC at ) depression (B f ) 



0.27 (<0.001 to 0.57) 0.001 (<0.001 to 0.24) 4.34 (3.78 to 4.70) 0.06 (<0.001 to 0.12) -2.41 (-5.55 to 1.07) 



all their known ancestors, comprised 2432 individuals. 
Mean pairwise k was 0.065 (median 0.061, inter-quartile 
range (IQR) 0.043-0.079) among all 2432 individuals 
and 0.069 (median 0.063, IQR 0.049-0.081) among 
the 2196 observed offspring (electronic supplementary 
material). Mean offspring / was 0.068 (median 0.064, 
IQR 0.038-0.090). 

A univariate animal model fitted to recruitment data for 
all 2196 offspring showed significant V A and h 2 (table 1). 
There was substantial inbreeding depression such that 
more inbred offspring were less likely to recruit, and 
males were more likely to recruit than females (table 1). 
Phenotypic recruitment probability did not differ between 
EPO and WPO across all 2196 offspring, and the EP 
status by sex interaction was not significant (table 1 and 
electronic supplementary material) . 

A bivariate animal model for male versus female recruit- 
ment showed that V A , h 2 and inbreeding depression in 
recruitment were similar in both sexes and that r mf w 1 
(table 2), indicating that genes underlying recruitment 
have similar additive effects in males and females. All off- 
spring were consequently pooled in subsequent bivariate 
analyses to maximize statistical power. 

(b) Male net paternity gain 

Male net paternity gain (AT E — N c ) was observed for 
293 individual males totalling 738 male-years during 



1993-2008. Of these males, 120 (41%) were observed 
in only 1 year (median 2 years, IQR 1-3, maximum 10). 
Net paternity gain varied from 11 to — 11, 6 to —6 and 
5 to — 5 through all offspring, sons and daughters respect- 
ively (electronic supplementary material) . The phenotypic 
correlation between net paternity gain through sons and 
daughters was 0.47 across all male-years. The pedigree, 
pruned to the 293 males and all their known ancestors, 
comprised 553 individuals. Mean pairwise k was 0.056 
(median 0.052 and IQR 0.023-0.077) among all 553 indi- 
viduals, and 0.071 (median 0.065 and IQR 0.048-0.085) 
among the 293 observed males. Mean male / was 0.057 
(median 0.047 and IQR 0.029-0.075). 

A univariate animal model fitted to male net paternity 
gain through all offspring indicated non-zero V A and h 2 , 
although the 95% CIs converged towards zero (table 3). 
Net paternity gain tended to decrease with increasing 
male /, but the 95% CI for inbreeding depression 
overlapped zero (table 3). 

(c) Genetic covariance 

The combined dataset comprising all adult males and off- 
spring observed during 1993-2008 comprised 2266 
individuals. The pedigree pruned to these 2266 individ- 
uals and all their known ancestors comprised 2445 
individuals. Mean pairwise k was 0.065 (median 0.061 
and IQR 0.043-0.079) among all 2445 individuals, and 
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0.068 (median 0.063 and IQR 0.048-0.081) among the 
2266 individuals with phenotypic data. 

Across all males and offspring, a bivariate animal model 
showed significant V A and h in both recruitment and male 
net paternity gain (table 4). Furthermore, there was sig- 
nificant negative genetic covariance between the two 
(table 4). Because £[iV E ] = 0.85, the mean estimated differ- 
ence in BV for recruitment between a female's EPO and the 
WPO they replaced was £[ABV] = -0.16 (95% CI -0.31 
to —0.01). Estimates of cov A (Ar E — N c ,W) were also nega- 
tive when analyses were repeated with net paternity gain 
through sons and daughters measured separately, although 
the 95% CI for sons overlapped zero (table 4). Because 
£[A?e S ] = 0.44 and £[iV EE) ] = 0.42, the mean estimated 
difference in BV for recruitment between a female's 
EP sons and the WP sons they replaced was 
£[ABV S ] = -0.15 (95% CI - 0.32 to 0.04), and the corre- 
sponding difference for a female's EP versus WP daughters 
was £[ABV D ] = -0.16 (95% CI -0.35 to -0.01). EPO 
therefore had a lower BV for recruitment than the WPO 
they replaced on average across all offspring, with similar 
estimated effects across sons and daughters. The estimate 
of £[ABV] w —0.16 equates to an average reduction in phe- 
notypic recruitment probability between maternal half-sib 
EPO and WPO of ca 0.02 in an average year. 

(d) Phenotypic variation 

A total of 1038 offspring hatched in 358 broods that con- 
tained at least 1 EPO at ringing, comprising 234 known 
mixed-paternity broods and 124 EPO-only broods. 
Across this restricted set of offspring, EP daughters 
tended to be less likely to recruit than WP daughters, 
while EP sons tended to be more likely to recruit than 
WP sons (figure 1). The EP status by sex interaction 
was significant and the estimated differences in recruit- 
ment probability were relatively large (figure 1). 
Although the 95% CIs for the sex-specific contrasts in 
recruitment probability between EPO and WPO over- 
lapped zero, they did not overlap the posterior modal 
difference for the opposite sex (figure 1). 

4. DISCUSSION 

One hypothesis explaining EPR by socially monogamous 
females is that a female's EPO are of higher additive gene- 
tic value than the WPO they replaced, and hence that 
EPR is under indirect selection through additive 'good 
genes' mechanisms [1-3]. This hypothesis differs from 
the hypotheses that EPR allows a female to produce 
additional offspring rather than solely replace WPO 
with fitter EPO (and is consequently under direct selec- 
tion), and that EPR is under indirect selection through 
'compatible genes' mechanisms [1,2,5]. However, the 
specific hypothesis that EPO are of higher additive gen- 
etic value than the WPO they replaced has not been 
explicitly tested in a wild population. The common 
short-cut approach of comparing phenotypes between 
observed maternal half-sib WPO and EPO from the 
same brood or litter is elegant in controlling for average 
phenotypic effects of environment and maternal genes 
but cannot distinguish paternal additive genetic effects 
from any individual environmental, parental or non-addi- 
tive genetic effects associated with paternity, and 
inevitably excludes EPO-only broods where EPR 
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Figure 1. Posterior modal estimates (and 95% credible inter- 
vals (CIs)) for the phenotypic probability of recruitment of 
extra-pair (EP) and within-pair (WP) daughters and sons 
estimated across 1038 known-sex offspring hatched in 
358 broods that contained at least 1 EPO. Posterior modes 
(and 95% CIs) for the main effects of inbreeding coefficient, 
EP status and EPO sex were -6.4 (-10.8 to -1.5), -0.58 
(-1.10 to -0.01) and -0.32 (-0.84 to 0.36), respectively. 
The EP status by sex interaction was significant (posterior 
mode: 0.86, 95% CI: 0.05-1.58). The posterior modes for 
the sex-specific contrasts in recruitment probability between 
EPO and WPO were -0.08 (95% CI -0.18 to 0.01) and 
0.05 ( — 0.04 to 0.15) for daughters and sons respectively. 



occurred but no WPO are available for phenotypic com- 
parison [2,8,10,13,14]. We suggest a direct approach to 
estimating the mean difference in additive genetic BV 
between EPO and the hypothetical maternal half-sib 
WPO they replaced (.EfABV]) by estimating the genetic 
covariance between fitness and male net paternity gain, 
and apply this approach to offspring recruitment in 
song sparrows. 

(a) Additive genetic (co)variances 
Because i:[ABV] is proportional to the genetic covariance 
between fitness and male net paternity gain, the existence 
of any indirect additive genetic benefit or cost of EPR 
through this route requires non-zero V A in both fitness 
and net paternity gain. Indeed, there was significant V A 
and h 2 in offspring survival to recruitment, one major fit- 
ness component. This was evident in both males and 
females, with an inter-sex genetic correlation of r mf w 1 . 
There was therefore some potential for female song spar- 
rows to alter the additive genetic value for recruitment of 
their offspring through social and/or EP mate choice. Fur- 
thermore, there was evidence of non-zero V A and h in 
male net paternity gain. Although 95% CIs converged 
towards zero in univariate analyses, estimates of V A and 
h 2 differed significantly from zero in bivariate analyses 
owing to increased power stemming from genetic covari- 
ance with recruitment. Indeed, the genetic covariance 
between net paternity gain and recruitment was signifi- 
cantly negative across all observed instances of EPR and 
across instances of EPR through daughters. The esti- 
mated effect was also negative across sons, although the 
95% CI overlapped zero. Opposite to widespread predic- 
tion, a female's EPO therefore had lower BV for 
recruitment than the WPO they replaced on average, 
representing an average additive genetic cost of EPR to 
polyandrous females. 



(b) Phenotype versus breeding value 

Across all 358 observed broods that contained at least 
1 EPO at ringing, EP daughters tended to be less likely 
to recruit than WP daughters, while EP sons tended to 
be more likely to recruit than WP sons with a significant 
EP status by sex interaction and substantial estimated 
effects (figure 1). Similar variation was observed across 
offspring hatched in known mixed-paternity broods 
[12]. EP sires therefore tended to produce daughters 
but not sons with low phenotypic recruitment probability. 
Such phenotypic analyses cannot determine whether 
observed variation reflects additive genetic, non-additive 
genetic and/or environmental differences between EP 
and WP sons and daughters; yet these different mechan- 
isms have very different implications for the evolution of 
female EPR [2,8,12]. Indeed, the estimated differences 
in mean BV for recruitment between EP and WP sons 
and daughters did not mirror observed phenotypic vari- 
ation. Observed EP sons tended to be more likely to 
recruit than observed WP sons (by ca 7% on average), 
even though EP sons tended to have lower BVs than the 
hypothetical WP sons they replaced (equating to a ca 
2% lower average genetic probability of recruitment). 
Observed EP daughters tended to be less likely to recruit 
than observed WP daughters (by ca 9%, on average), 
while EP daughters had a ca 2% lower average genetic 
recruitment probability than the hypothetical WP daugh- 
ters they replaced. These discrepancies between estimated 
phenotypic and additive genetic effects should not be 
over-interpreted because BVs but not phenotypes can be 
measured for hypothetical WPO that did not exist, and 
the 95% CIs for the estimated sex-specific phenotypic 
effects overlapped the estimated genetic effects. However, 
the implication is that there may be sex-specific environ- 
mental, parental and/or non-additive genetic influences 
on the relative recruitment of EPO versus WPO [12]. 
One conceivable mechanism could be that EPO tend to 
hatch earlier within a brood than WPO [8,14] and early 
hatching benefits sons but not daughters. Female EPR, 
and hence phenotypic differences between EPO and 
WPO, could also arise because EP sires are less closely 
related to or otherwise more 'compatible' with a female 
than her socially paired male [1,2]. Substantial inbreeding 
depression was evident in offspring recruitment in both 
sexes, implying non-additive genetic variation (directional 
dominance) in recruitment (electronic supplementary 
material). A further conceivable mechanism could there- 
fore be that females avoid inbreeding through EP sons 
more than daughters. These hypotheses remain to be 
tested. However, irrespective of the mechanisms, our ana- 
lyses demonstrate that phenotypic comparisons among 
observed maternal half-sib EPO and WPO may not 
exactly mirror the sex-specific difference in BV between 
EPO and the WPO they replaced, either quantitatively 
or qualitatively. Such phenotypic comparisons should 
therefore be interpreted with some caution in the context 
of inferring indirect selection on female EPR through 
increased additive genetic value of EPO compared with 
their WPO maternal half-sibs [3,8,13,14]. 



(c) Evolution of female extra-pair reproduction 

A model predicting the rate of evolution of female 
EPR due to indirect selection through increased 
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additive genetic value of EPO has been derived as 
= ^pEPO ' °pEPO " 4w, where Aj is the number of phe- 
notypic standard deviations by which the mean EPR 
rate would evolve in one generation through such indirect 
selection alone, ^p EPO is the heritability of the proportion 
of a female's offspring that is sired by an EP male, cr pEPO 
is the phenotypic standard deviation of this proportion 
and <i EW (strictly defined) is the difference in paternal 
additive genetic value for fitness across all EPO and 
WPO [3]. Values of h 2 pEPO « 0.12 and cr pEPO « 0.32 
have been estimated for Mandarte's song sparrows [27] . 
The value of d EW can be approximated from field studies 
as the mean difference in phenotypic fitness between 
maternal half-sib EPO and WPO [3]. However, E[AEV] 
may provide a better approximation because it eliminates 
any differential parental or environmental effects on phe- 
notype. Taking £[ABV] for recruitment w —0.16 gives an 
approximate point estimate of Ai « —0.006. More com- 
prehensive estimation of Ai using this approximation 
would require £;[ABV] to be estimated for total fitness 
from conception rather than solely recruitment. This is 
challenging because the full distribution of fitness is diffi- 
cult to measure and often highly skewed, and r mf for total 
fitness may not be one. However, recruitment is a major 
determinant of fitness in song sparrows and other ver- 
tebrates [20,32]. Therefore, unless EPO have much 
higher additive genetic value for survival to hatch and/or 
reproductive success than the WPO they replaced, our cur- 
rent data indicate that A : may be slightly negative in song 
sparrows, implying evolution towards less EPR. Other 
explanations for female EPR, including direct selection 
and other components of indirect selection, may therefore 
be required [2-6,26] . The observation that male net pater- 
nity gain (or loss) through EPR is significantly heritable is 
interesting in this context; genetic covariances between 
net paternity gain and other components of male and 
female fitness now need to be quantified. 

(d) The value of maternal half-sib comparisons 

Comparing BVs between EPO and the hypothetical 
maternal half-sib WPO they replaced is insightful in 
explicitly indicating the additive genetic benefit that polyan- 
drous females gain from EPR, including for EPO-only 
broods, and in considering the degree to which standard 
phenotypic comparisons between maternal half-sibs may 
accurately indicate this benefit. However, ultimately, the 
hypothesis that female EPR is under indirect selection 
through additive genetic mechanisms predicts and requi- 
res covariance between genes underlying female EPR and 
paternal and/or offspring fitness arising from within-pair 
and extra-pair reproduction across all polyandrous and 
monogamous females [3,9] . Estimation of this covariance 
has been simplified to a (phenotypic) comparison between 
half-sib EPO and WPO of polyandrous females purely for 
empirical tractability in wild populations [3, 1 0] . However, 
EPO typically occur in a minority of broods in socially mon- 
ogamous birds ([2]; ca 44% in Mandarte's song sparrows 
[21]). Even explicit comparisons of BV between maternal 
half-sib EPO and WPO may poorly indicate the overall gen- 
etic covariance between female EPR and offspring or 
paternal fitness if WPO-only broods occur non-randomly 
with respect to genetic values for female EPR, offspring 
and/or social male fitness (as might be hypothesized). 



WPO-only broods cannot contribute to phenotypic or 
quantitative genetic comparisons between observed or 
hypothetical maternal half-sib EPO and WPO because no 
EP sire exists or is known. Field studies aiming to dissect 
the selective forces driving female EPR, and definitively esti- 
mate <i EW and Aj, must therefore go beyond the maternal 
half-sib comparisons that have become the default approach 
and estimate overall genetic covariances among EPR and 
paternal and offspring fitness. 

We thank everyone who contributed to long-term data 
collection on Mandarte, Peter Arcese, Lukas Keller and 
Erik Postma for valuable comments, Jarrod Hadfield for his 
generous insights and assistance with deriving the 
expression for .E[ABV], and the UK Royal Society and 
NERC for funding. 



REFERENCES 

1 Jennions, M. D. & Petrie, M. 2000 Why do females mate 
multiply? A review of the genetic benefits. Biol. Rev. 75, 
21-64. (doi:10.1017/S0006323199005423) 

2 Griffith, S. C, Owens, I. P. F. & Thuman, K. A. 2002 
Extra pair paternity in birds: a review of interspecific vari- 
ation and adaptive function. Mol. Ecol. 11, 2195-2212. 
(doi: 10. 1046/j.l365-294X.2002. 01613.x) 

3 Arnqvist, G. & Kirkpatrick, M. 2005 The evolution of 
infidelity in socially monogamous passerines: the strength 
of direct and indirect selection on extrapair copulation 
behavior in females. Am. Nat. 165, S26-S37. (doi: 10. 
1086/429350) 

4 Akcay, E. & Roughgarden, J. 2007 Extra-pair paternity in 
birds: review of the genetic benefits. Evol. Ecol. Res. 9, 
855-868. 

5 Westneat, D. F. & Stewart, I. R. K. 2003 Extra-pair pater- 
nity in birds: causes, correlates and conflict. Annu. Rev. 
Ecol. Evol. Syst. 34, 365-396. (doi: 10.1 146/annurev. 
ecolsys.34.01 1802.132439) 

6 Griffith, S. C. 2007 The evolution of infidelity in socially 
monogamous passerines: neglected components of direct 
and indirect selection. Am. Nat. 169, 274-281. (doi: 10. 
1086/510601) 

7 Griffith, S. C. & Immler, S. 2009 Female infidelity and 
genetic compatibility in birds: the role of the genetically 
loaded raffle in understanding the function of extra-pair 
paternity. J. Avian Biol. 40, 97-101. (doi: 1 0. 1 1 1 
1600-048X.2009. 04562.x) 

8 Magrath, M. J. L., Vedder, O., van der Velde, M. & 
Komdeur, J. 2009 Maternal effects contribute to the 
superior performance of extra-pair offspring. Curr. Biol. 
19, 792-797. (doi:10.1016/j.cub.2009.03.068) 

9 Kirkpatrick, M. & Barton, N. H. 1997 The strength of 
indirect selection on female mating preferences. Proc. 
Natl Acad. Sci. USA 94, 1282-1286. (doi:10.1073/ 
pnas. 94. 4.1282) 

10 Sheldon, B. C, Merila, J., Qvarnstrom, A., Gustafsson, 
L. & Ellegren, H. 1997 Paternal genetic contribution of 
offspring condition predicted by size of male secondary 
sexual character. Proc. R. Soc. Lond. B 264, 297-302. 
(doi: 10.1 098/rspb. 1 997.0042) 

1 1 Schmoll, T, Schurr, F. M., Winkel, W., Epplen, J. T. & 
Lubjuhn, T. 2009 Lifespan, lifetime reproductive 
performance and paternity loss of within-pair and extra- 
pair offspring in the coal tit Periparus ater. Proc. R. 
Soc. B 276, 337-345. (doi:10.1098/rspb.2008.1 1 16) 

12 Sardell, R. J., Keller, L. F, Arcese, P. & Reid, J. M. 
2011 Sex-specific differential survival of within-pair and 
extra-pair offspring in song sparrows. Proc. R. Soc. B 
278, 3251-3259. (doi:10.1098/rspb.201 1 .0173) 



Proc. R. Soc. B (2012) 



1708 J. M. Reid & R. J. Sardell Genetic value of extra-pair offspring 



13 Johnson, L. S., Brubaker, J. L., Johnson, B. G. P. & 
Masters, B. S. 2009 Evidence for a maternal effect bene- 
fiting extra-pair offspring in a songbird, the house wren 
Troglodytes aedon. J. Avian Biol. 40, 248-253. (doi:10. 
1 1 1 l/j.l600-048X.2009. 04777.x) 

14 Ferree, E. D., Dickinson, J., Rendell, W., Stern, C. & 
Porter, S. 2010 Hatching order explains an extrapair 
chick advantage in western bluebirds. Behav. Ecol. 21, 
802-807. (doi:10.1093/beheco/arq056) 

15 Lynch, M. & Walsh, B. 1998 Genetics and analysis of 
quantitative traits. Sunderland, MA: Sinauer. 

16 Kruuk, L. E. B. 2004 Estimating genetic parameters in 
natural populations using the 'animal model'. Phil. 
Trans. R. Soc. Lond. B 359, 873-890. (doi:10.1098/ 
rstb.2003.1437) 

17 Postma, E. 2006 Implications of the difference between 
true and predicted breeding values for the study of 
natural selection and micro-evolution. J. Evol. Biol. 19, 
309-320. (doi:10.1111/j.l420-9101. 2005. 01007.x) 

18 Hadfield, J. D., Wilson, A. J., Garant, D., Sheldon, B. C. & 
Kruuk, L. E. B. 2010 The misuse of BLUP in ecology 
and evolution. Am. Nat. 175, 116-125. (doi:10.1086/ 
648604) 

19 Kruuk, L. E. B. & Hadfield, J. D. 2007 How to separate 
genetic and environmental causes of similarity between 
relatives. J. Evol. Biol. 20, 1890-1903. (doi: 1 0. 1 1 1 1/j. 
1420-9101. 2007. 01377.x) 

20 Smith, J. N. M., Keller, L. E, Marr, A. B. & Arcese, P. 
2006 Conservation and biology of small populations: the 
song sparrows of Mandarte Island. New York, NY: 
Oxford University Press. 

21 Sardell, R. J., Keller, L. E, Arcese, P., Bucher, T. & Reid, 
J. M. 2010 Comprehensive paternity assignment: geno- 
type, spatial location and social status in song sparrows, 
Melospiza melodia. Mol. Ecol. 19, 4352-4364. (doi: 10. 
1 1 1 1/j. 1365-294X.2010. 04805.x) 

22 Keller, L. E, Jeffery, K. J., Arcese, P., Beaumont, M. A., 
Hochachka, W. M., Smith, J. N. M. & Bruford, M. W. 
2001 Immigration and the ephemerality of a natural 
population bottleneck: evidence from molecular markers. 



Proc. R. Soc. Lond. B 268, 1387-1394. (doi:10.1098/ 
rspb.2001.1607) 

23 Postma, E., Heinrich, E, Roller, U., Sardell, R. J., Reid, J. 
M., Arcese, P. & Keller, L. F. 20 1 1 Disentangling the effect 
of genes, the environment and chance on sex ratio variation 
in a wild bird population. Proc. R. Soc. B 278, 2996-3002. 
(doi:10.1098/rspb.2010.2763) 

24 Keller, L. F. 1998 Inbreeding and its fitness effects in an 
insular population of song sparrows (Melospiza melodia). 
Evolution 52, 240-250. (doi:10. 2307/2410939) 

25 Reid, J. M. & Keller, L. F. 2010 Correlated inbreeding 
among relatives: occurrence, magnitude and impli- 
cations. Evolution 64, 973-985. (doi: 1 0. 1 1 1 1/j. 1 558- 
5646.2009.00865.x) 

26 Reid, J. M., Arcese, P., Sardell, R. J. & Keller, L. F. 201 1 
Additive genetic variance, heritability and inbreeding 
depression in male extra-pair reproductive success. Am. 
Nat. Ill, 177-187. (doi:10. 1086/657977) 

27 Reid, J. M., Arcese, P., Sardell, R. J. & Keller, L. F. 201 1 
Heritability of female extra-pair paternity rate in song 
sparrows (Melospiza melodia). Proc. R. Soc. B 278, 
1 1 14- 1 120. (doi: 10. 1098/rspb.2010. 1704) 

28 Reid, J. M., Arcese, P. & Keller, L. F. 2006 Intrinsic 
parent -offspring correlation in inbreeding level in a 
song sparrow (Melospiza melodia) population open to 
immigration. Am. Nat. 168, 1-13. (doi:10. 1086/504852) 

29 R Development Core Team 2009 R: a language and 
environment for statistical computing. Vienna, Austria: R 
Foundation for Statistical Computing. 

30 Hadfield, J. D. 2010 MCMC methods for multi-response 
generalized linear mixed models: the MCMCglmm R 
package. J. Stat. Soft 33, 1-22. 

31 Nakagawa, S. & Schielzeth, H. 2010 Repeatability for 
Gaussian and non-Gaussian data: a practical guide for 
biologists. Biol. Rev. 85, 935-956. (doi: 1 0. 1 1 1 1/j. 1469- 
185X.2010. 00141.x) 

32 Gaillard, J. M., Festa-Bianchet, M. & Yoccoz, N. G. 1998 
Population dynamics of large herbivores: variable recruit- 
ment with constant adult survival. Trends Ecol. Evol. 13, 
58-63. (doi:10.1016/S0169-5347(97)01237-8) 



Proc. R. Soc. B (2012) 



